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The computer program 'Fecom.nb' implementing the Fedosov *-product in Darboux coordinates 
is presented. It has been written in Mathematica 6.0 but it can be easily modified to be run in some 
earlier version of Mathematica. To optimize computations elements of the Weyl algebra are treated 
as polynomials. Several procedures which order the terms are included. 



^! 1 Introduction 

» 

' A phase space formulation of quantum mechanics was introduced by Moyal [1] in the middle of the 20th 

century. Moyal's way of thinking about quantum world arose from earlier works by Weyl [5], Wigner 
[3] and Groenewold [4j. The phase space quantum mechanics in its original version was restricted to 
I systems with the phase spaces M^". It was a complete physical theory i.e. formulas defining a *-product 

^ ' of functions, mean values, eigenstates and describing a time evolution were given. 

As far as we know an extension of the Moyal quantum mechanics on nontrivial phase spaces containing 
a quantum state representation and a method of searching for eigenvalues and eigenstates does not exist. 
^f) • The starting point for constructing a theory of this kind must be some modification of the Moyal *- 

product. 

, The concept of a *- product on an arbitrary phase space was proposed by Bayen ei a/ [5], [6j. In 

this approach quantum theory was considered as a deformation of classical mechanics with the Planck 
constant ft, as a deformation parameter. 

According to Bayen et al the ^-product of smooth functions A and B defined on a phase space can 
J\. • be written as follows 

?-H ' X +°° 

_Cd ; A*B = A-B+—{A,B}p + J2^''^k{A,B), (1.1) 

fc=2 

where A is a deformation parameter, {•,•}_? denotes the Poisson bracket and -Dfe's, 2 < k are some 
bidifferential operators. Works by Bayen et al established only general frames of deformation quantization 
and did not contain explicit formulas describing operators Dk so each case required individual analysis. 

In 90's of the previous century two outstanding papers devoted to a practical realisation of the Bayen 
et al quantization program on a symplectic manifold endowed with a symplectic torsionfree connection 
were published by Fedosov [7j , [8j . He found a purely algebraic construction of the ^-product determined 
by this symplectic connection. Since every symplectic manifold may be equipped with a symplectic 
connection, his method can be applied on any symplectic manifold. 

Explicit form of the *-product in this formalism is extremely complicated. However the Fedosov 
algorithm is based on two recurrence relations and therefore the ^-product of two functions A and B 
can be found up to an arbitrary fixed power of the deformation parameter after a finite number of steps. 
Moreover, his H<-product can be computed by a computer program. In this paper we present a program 
written in Mathematica 6. To improve efficiency of the code we restrict ourselves to Darboux coordinates. 
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Our contribution is a 'kernel' on which procedures for individual purposes can be built. In the 2nd 
section of this paper we show, how elements of the Fedosov quantization program have been implemented 
in the program. The next section is a short instruction for users. Some tests end the paper. 

In all formulas, in which summation limits arc obvious, we use the Einstein summation convention. 

Names of procedures are written in a bold typestyle. Names of variables are in italic. 



2 The Fedosov's construction of deformation quantization 

At the beginning we explain how the Fedosov algorithm has been implemented in our computer program. 
We assume that the Reader is familiar with the Fedosov method. For details see [8]. 

Let {M,uj,r) be a 2ti-D Fedosov manifold and A ~ {{Uz,4>z)}zeJ a Darboux atlas on this manifold. 
By to we mean the symplectic 2-form. In each chart belonging to the atlas A 

n 

CO = "^dx' Adx'+''. (2.1) 
i=i 

The symbol F denotes a symplectic connection i.e. torsionfree connection preserving the symplectic 2- 
form. Locally in Darboux coordinates the symplectic connection is determined by a set of coefficients 
{^ijk "== r'^w,;;} symmetric in all indices [71 [H [9]. 

Let h denote a deformation parameter. We assume that it is positive. By y^, . . . we denote the 
components of an arbitrary vector y belonging to the tangent space T^Ai at the point -p G Ai with respect 
to the natural basis (^)p determined by the chart {Uz,4>z) such that p GHz- 

We introduce the formal series 



oo 

dcf. 

a ~ 



Y^h^a^,,,...,,V^' ...V^\ k>0 (2.2) 

/=o 

at p. For / = we put a = h^au- By ak.ji...ji we denote the components of a covariant tensor symmetric 
with respect to the indices {ji , . . . ,ji} taken in the basis dx^^ ... da;-" . 

In our program the deformation parameter h is denoted by h. The components , . . . , are elements 
of a vector Yy, so in the computer program is represented by Yy[[i]]=y[i]. 

The part of the series a standing at h'' and containing I components of the vector y will be denoted 
by a[k,l]. Thus 

oo oo 

a = ^^/i'=a[fc,/]. (2.3) 

fe=0 1=0 

The degree deg(a[fc, I]) of the component a[k, I] is equal to 2k + I. 

Notice that since akj-^...j, are completely symmetric in the indices {ji, . . . , j;}, the element a defined 
by the formula (|2.2p can be understood as the polynomial 



« = E E (y'r ■ ■ ■ (y'")^^" ' (2-4) 

z=0 fe=0 

where 

< ii, . . . , i2n < z — 2k , ii + ■ ■ ■ + i2n = z — 2k. 

The symbol [|] denotes the floor of |. The relation between the tensor components ak,ji...ji and the 
polynomial coefficients dk,ii...i2n reads 

(2.5) 
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Using the polynomials instead of the completely symmetric tensors shortens drastically the time required 
for computations. 

Let Pp*7W [[ft,]] denote a set of all elements a of the form (|2.2p at the point p. The product o : Pp* [[/i]] x 
P*M[[h]] P*M[[h]] of two elements a,b e P*M[[h]] is the mapping 

aob = } -r -TT • • • w 7-- TT- TT—, (2-6) 

^ i! V 2 y . . . dy'* dy^^ . . . dy->* ^ ' 

where u)^'^ is the inverse tensor to ujkj i-e. 

The pair (Pp*A^[[/i]], o) is a noncommutative associative algebra called the Weyl algebra. All a,b G 
(Pp*A^[[ft]],o) fulfill the relation 

deg(ao6) deg(a) +deg(6). (2.7) 

Our program works in a Darboux chart. In this chart formulas defining the o-product are simple and 
we may treat elements of the Weyl algebra as the polynomials instead of the tensors. Moreover, the 
symplectic connection is represented by an element of the Weyl algebra (see a relation (|2.16p ) and we 
calculate exterior covariant derivatives in terms of the Weyl algebra (|2.15|1 . The o-product is computed 
in three steps. 

1. The module Circl finds the o-product of the monomials [y^Y {y'^^^Y o {y'^Y {y'^^^Y , 
1 < i < n, < r,j,s, k. As it has been proved [TUj 

min[r,fc]+min[j,s] / .,\ t 

(y'r(2/'+")^°(y')^(2/'+")'=Hj!s!fc! Yl V (2/'r-'^"*(2/'+")'-*"'"*x 

t=o \ / 

min[j,s,i] ^ 

' + a)! ij ^ ay. is- ay. ik^t + «)! " ^'''^ 

2. The o-product of the monomials 

{{y'YHy'r ■ ■ ■ iy'"Y'") ° {{y^YHyy^ • ■ • iy'"Y'") - 
{{y'r (y»+i)«"+i o {yy^ (2/"+!)^"+^) ■ • • [{y-Y" (y'")'=" ° ify- (2/'"^") (2.9) 

is generated in the procedure Circ2. 

3. Using bilinearity of the o-product the module Circproduct finds the product of the two polynomials 
from the Weyl algebra (Pp*A^[[ft]], o). 

In the operation of the o-multiplication of the polynomials aob most of time is used to compute the 
product (j2.9p of the monomials. Hence we prepared the procedure Ord, which orders the polynomial in 
y^, . . . , y^" in such a way that every term [y^Y^ {y'^Y^ ■ ■ ■ (y^")*^" appears only once. For example 



Ord 



1 1 



ay[l] + f[x[l],x[2]]y[l] + -y[2] + -y[l] = - + a + /[x[l],x[2]] h/[l] 



2"' ' 3 
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Taking a set of the Weyl algebras (Pp*A^[[/i]], o) at all points of the manifold M we constitute the 
Weyl bundle 

V*M[[h]]^^- U {P;M[[h]],o). 
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Geometrical structure of the Fedosov deformation quantization is based on an m-differential form 
calculus with values in the Weyl bundle. Locally such a form can be written as follows 



h''ak,j,...j,,s^...sAx\ . . . . ..y^^dx'^ A • • • A dx'"^ , (2.10) 

1=0 

where < m < 2n. Now ak.ii...ii,ji...j^(x^ , ■ ■ ■ ,2^^") are components of smooth tensor fields on JV[ and 

C°°(TA^) 3 y '°="^ is ^ smooth vector field on AA. We use the same symbol for the vector field 

y G C°°{TM) and the vector y G TpA^. In the program we treat these objects as entries of the same 
vector Yy. 

From now on we will omit the variables [x^ , . . . , a;^") in ak,ji...ji.si...s^{x^ , ■ • ■ , a;^"). 

The forms dUHl) are smooth sections of the direct sum V*M[[h]] ® A ®^=o(^*-^[[^]] ® A™). 
By A™ we mean a smooth field of m- forms on the manifold M. The o-product a o 6 of two forms 
a e C°°{V*M[[h]] (g> A"!) and b G C'^{r*M[[h]] ® A™^) is understood as the o-product of the Weyl 
algebra elements and the exterior product 'A' of the forms. 

Computations in the Fedosov deformation quantization are done on 0-forms and 1-forms. Hence we do 
not prepare a general module calculating the exterior product of forms. We give only the A-multiplication 
of 1-forms dx' A dx^ , 1 < i,j < 2ti. 

In our program the elements of a local basis of A^ are entries of a vector dxx written as dx[i]. The 
product dx* A dx^ has been denoted as dxx[[i]]**dxx[[j]]. The symbol '**' is used in Mathematica for 
some associative but noncommutative product. 

For our purposes it is sufficient to implement antisymmetry of the A-product. Hence we declare values 
of the products dx' A dx^ for I < i < j < 2n first and then introduce antisymmetry. The exterior product 
dxx[[i]] A dxx[[j]] is denoted as Wedgeproduct[dxx[[i]], dxx[[j]] ], where 

Vi<i<j<2ji Wedgepr o duct [ dxx [ [i] ], dxx [[j]] ] dxx[[i]] * *dxx[[j]] 

and 

Wedgeproduct[xl9_,xl9 ] := , Wedgeproduct[xl3 ,yl3 ] := Wedgepr o duct [y 13, x 13]. 
Therefore for example 

Wedgeproduct[dxx[[l]],dxx[[2]] ] = dx[l] * H<dx[2], 

Wedgeproduct[dxx[[l]],dxx[[l]]] = 0, 

Wedgeproduct[dxx[[2]],dxx[[l]] ] = -dx[l] * *dx[2]. 

The o-product of 1-forms a G V*M[[h]] ® A^ and b G ® A^ with values in the Weyl algebra is 

computed by the module Fullproduct. 

An important role in our considerations is played by a commutator of forms. By definition the 
commutator of forms a G C°° {V* M[[h]] A™i) and b G {V* M[[h]] A™^) is the form [a, 6] G 

C'^{V*M[[h]] ® A"'i+™^) 

[a,b] '^^^ ao6- (-l)'"i-"^6oa. (2.11) 

In the program the commutator of 1-forms a and b is calculated by the module Commutatorll. 
The operator S'^ : C°°{V*M[[h]] ® A") ^ C°°{r*M[[h]] ® A™-i) is defined by 

TT— W*" TTT I a for / + ?71>0, 

6-^a={ (2.12) 
for / + m = 0, 

where I is the degree of a in y^ 's i.e. the number of y-^'s. The operator raises the degree of the forms 
of P*M[[h]]A in the Weyl algebra by 1. 
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The operator 5 ^ is linear so we need to implement its action on the monomials (y^)*^ . . . (y^")'^"(ia;-' 
and (y^y^ ■ ■ . {y'^")''^" dx^ A dx'^ . By direct calculations it can be checked that 

S-'Uy^^ . . . {y'^T"dxA = — ]. ^ J iv'T' ■ ■ ■ {y'n'^"y') (2.13) 

V /?! + ...+ l2n + 1 ^ ^ 

and 

d-U{y'r • ■ • {y^n^-^-dx^ Adx^) =— ^ ^o i^y'y" ■ ■ ■ {y'n'^-y'dx^ - [y'r ■ ■ ■ {y^n'^-y^dxA. 

V /«! + ...+ l2n + 2 V / 

(2.14) 

Acting of the operator S ^ on 2-forms is realized by the module Deltal and on 1-forms by DeltaO. 
Before using these modules the forms are ordered by Ordnung2 and Ordnungl respectively so each 
element {y^Y^ . . . (y'^^Y^" dx^ A dx'' or {y^Y^ ■ ■ ■ {y'^")''^" dx^ appears only once. 

The exterior covariant derivative of a form a G C°°('P*Al[[/i]] A™) determined by a connection 
l-form7 G C'^{P*M[[h]]^A^) is the linear operator 9^ : C°°(P*X[[/i]] A™) -> C°°{P*M[[h]]^A'^+^) 
defined in a Darboux chart by the formula 

dja'^=' da+ ^[jja]. (2.15) 



In the case when 7 represents the induced symplectic connection we put 

7'™=°'r = ir,,,2/Vd^^ (2.16) 



The 1-form of the symplectic connection F is called symplecticconnectionlform and it is computed by a 
command of the same name. 

The curvature form of a connection 1-form 7 in a Darboux chart can be expressed by the formula 



z 1 

i?^ = ^7+ — [7,7] = d7+ -707. (2.17) 



In the program a curvature of some connection kona is computed by the function Curvaturecalcula- 
tion. Hence the symplectic curvature 2-form is found as Curvaturecalculation[symplecticconnectionlform] . 

The crucial role in the Fedosov deformation quantization is played by an Abelian connection F. The 
Abelian connection proposed by Fedosov is of the form 

r ^ ujijy'dx^ + T + r, (2.18) 

where the series r is at least of the degree 3 and it is defined by the recurrence relation 



S-'Rr + 5-' [drr + -ror] . (2.19) 



Let r[z] denote the component r[z] X]fc=o ' h^'''''m[2k, z — 4k]dx"^, z > 3 of r of the degree z. As it 
was shown [11], [12] 

r[3] = 6-'R.T, 

r[z] = S-'^ ^rr[z - 1] + o r[z + 1 - j]j , z > 4. (2.20) 

The Abelian connection is written as the vector matrixAbelconnection. The element of the degree z 
is saved as matrixAbelconnection[fz]]. There is only one exception. Since the term ujijy^dx^ does not 
appear in further operations, we put matrixAbelconnection[[l]] = 0. The quantities matrixAbelconnec- 
tion[[z]] for 3 < z are calculated by the function CorrAbelconnection and ordered by the procedure 
Ordnungl before saving. The Abelian connection in the form of series is generated by the procedure 
Abelianconnectionseries. 
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The subalgcbra r*M[[h]]f, C C°° {V* M[[h]] (g> A") consists of flat sections i.e. such that df,a = 0. For 
any qq € C°°{Ai) there exists a unique section a e 7'*A^[[ft,]]j=, such that the projection 

, ^ dcf. I 

a[a) = a|y=o = oo- 
The element a = a^^{ao) can be found by the iteration 



ao + S^^ \^dTa+ j^[r,a]j . (2.21) 
This relation means that 

a[0] = ao, 

. z-2 

a[z]^6-^(dra[z-l] + ^^\^r[z + l-l],a[l]y z>l. (2.22) 

1=1 

Remember that r[l] = r[2] = 0. The series representing multiplied functions A and B are written as the 
vectors matrixA, matrixB respectively. The part of degree z is on the (z + l)'th position of the vector. 
The term a[z] is generated by the module Flatcorrection and it is ordered by the function Ord. 

Using the one-to-one correspondence between 'P*A^[[/i]]p and C°°{A4) we introduce the associative 
star product of functions ao,bo E C^{M) 

ao * bo fT ((T-i (ao) oa-\bo)). (2.23) 

The *- product ()2.23|) with A = —ft, is of the form and in physics it is interpreted as the quantum 
multiplication of obscrvablcs. Notice that in the Fcdosov notation the correspondence principle reads 



{ao,bo}p ^J^{°'° * bo - bo* ao^ 



There is no need to calculate the complete product a ^(aq) o a ^{bo) and then to project it on 7W. 
The part of the product ao * bo standing at h'^ must be the projection of X]i^=o '''"^ ('^o) ['] ° <J^^{bo)[2k — I] . 
Since we are looking for the term at h'^ we consider only elements of this form. Moreover, for / > each 
expression a~^{ao)[l] or a^^(6o)[2A; — I] contains at least one y so we deduce that 



ao*bo^ao-bo + J2^[ J2 ^^-oM ° cj-\bo)[2k \ . 

k=l \ 1 = 1 ) 

Among the products (T~^(ao)[^] o CT~^(6o)[2fc — l\ only these which do not generate any y's are important. 
Therefore, for a fixed /, we check which of the series: tT~^(ao)[Z] or cr~^(&o)[2fc — /] is shorter. Then for 
each element of this shorter series we look for its counterpart giving in the o-product a term containing 
h!^ and then multiply these elements. 

Let the term 

A = h^- A{x{\l . . . , x{2n\) {y^r- [y^f^ . . . (j/^")^- 
from (7^^(ao)[2s + ii + . . . + i2„] be chosen. Its counterpart 

B = ft^— . B{x[ll . . . , x[2n]) {y'Y-^^ {y^f-^^ . . . (y^")*" 

belongs to a^^{bo)[2k — 2s — ii — ... — i2„]. Then using (|2.8p and (|2.9p one gets 

a{A oB)^h''- Aix[l],. . . , x[2n]) ■ B(a-[1], . . . , x[2n])(-l)'-+i+-+'-" ^ h\ ■ ■ ■ . 

The final product is represented by the vector productmatrix. The terms containing are on the 
(fc-l-l)'th position. These terms are computed by the module Rdegreecorrection. Appriopriate elements 
are found by the function Lookingfor and the reduced o-product is calculated by Partialcirc. 
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3 The program description 



After starting the program the User is asked to give n i. e. one half of the dimension of the phase 
space {A4,uj). The number n is usually the dimension of the configuration space. Then the program 
reads coefficients of the symplectic connection F. The symbol Gamma[i,j,k] denotes the coefficient Tijk- 
Since Tijk arc symmetric in all indices, the computer asks only for the components Tijk such that 
^<'i^j<k< 2n. In the program local coordinates on A4 are elements of the vector xx. The i-th 
coordinate is then xx[[i]] but, since by the definition xx — Array [x, 2n], the shorter form x[i], i = 1, . . . ,2ti 
is acceptable. 

In the next step it is required to declare the range of approximation of the *- product. In the case 
when the multiplied functions do not contain the deformation parameter h, the range is just a maximal 
power of h appearing in the product. This quantity, denoted by hpower, determines the number of 
repetitions in the recurrence procedures ([2?20l) and ([2?22l) . Indeed, for a given hpower it is necessary to 
make (2 • hpower — 1) repetitions. 

Finally we declare the functions A and B to be multiplied. 

To enable the User following the Fedosov algorithm, the program prints in an expanded form the I- 
form of the symplectic connection F, its curvature i?r, the sum of the symplectic connection and the 
Abelian correction r and the flat sections a-~-^{A) and a~^{B). All these quantities are known to the 
degree (2 • hpower — 1). 

The program do not check if the components of the symplectic connection Tijk and the multiplied 
functions A and B are differentiablc sufficiently many times. Please, be careful and do not use symbols 
which may be confused with Mathcmatica commands. 

Some basic information about variables and modules can be seen after the use of command ?name- 
ofproccdure. A table below contains a list of symbols used in our paper and their counterparts in the 
program. 



Article 


program 


h 


h 




x[i] or xx[[i]] 


y' 


y[i] or Yy[[i]J 




dx[i] or dxx[[i]] 
dxfi] ** dx[j] or dxx[[i]] ** dxx[[j]] 


dx'- A dx^ 


Tijk 


Gamma[i,j.k] and symplecticconnection[[i, j, k]] 


T 


symplecticconnectionlform 


Rr 


symplecticcurvature 


r[z], z > 3 


matrixAbelconnection[fz]] 


a[z], 6[z], 2 > 


matrixA[[z+l]], matrixBf[z+l]] 



4 Tests 

To check efficiency of our program we applied it to calculate three examples of the Fedosov ^-product. 
The computations were done in Mathematica 6.0.1.0 on the platform Microsoft Windows XP Professional 
x64 version 2003. The computer used the processor AMD Athlon 64x2 Dual 4400+ and it was equipped 
with 2GB of RAM. During the tests Mathematica was the only one running program. 

1. The Moyal product for rt = 1. 

This is a basic example of the Fedosov *- product on the flat manifold M^. In the Cartesian co- 
ordinates the symplectic connection 1-form F = 0. Hence the symplectic curvature vanishes and 
the same happens to the Abelian correction r. We multiplied two functions ^(^^[[l]], a;a;[[2]]) * 
z)(a;a;[[l]], xa;[[2]]) up to the terms standing at h'^. Using the function Timing we saw that the time 
used for generation the symplectic connection 1-form, its curvature and the Abelian correction was 
negligible. The matrix representing the flat section a^^ ( i(;(xa;[[l]] , a::x[[2]]) j was built in 0.219 s and 
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the matrix containing a ^ yv{xx\\l\\,xx[[2]])j in 0.203 s. The ^-product was found in 0.016 s. 

2. The general form of the Fedosov ^-product on a 4- dimensional Fedosov space. 
The expression representing a symplectic connection 1-form is 

4 

r- 51 9^3k{x[\],x[2],x[?>\,x[^)y\l]y[i\dx[k]. 

i,j,k=l 

The symplectic connection 1-form was formed by 0.016 s. Its curvature 2-form was computed in 
0.484 s. We searched the *-product up to of two functions w(a;[l], a:[2], a;[3], a;[4])*?;(a;[l], a;[2], a;[3], a;[4]) 

The correction r was built in 14 s. The series ^ui(x[l], a;[2], a;[3], a;[4])^ and ^w(a;[l], a;[2], a;[3], a::[4]) 

in an expanded form were generated by approx. 180 s each. Their ^-product was found in 0.046 s. 

3. The symplectic space (R^,w) with the symplectic connection determined by coefficients 

Fin = -x[2] , Fna = , Fiaa = , F222 = 
in a global Darboux chart. 

We generated the eigenvalue equation for a: [2] up to h^. This problem was solved in terms of the 
Fedosov formalism in |13j . Opposite signs of Fm in [13j and in our paper were the consequence of 
different sign conventions. The Wigner eigenfunction of x[2] was denoted by w[a;[l], a;[2]]. 

The symplectic connection 1-form F = — ix[2]y[l]^c?.T[l] and the symplectic curvature 2-form 
Rr = ^j/[l]^da;[l] A dx[2]. The time used for finding these two expressions was neghgible. The 
correction r to the Abelian connection was calculated in 0.016s, the flat sections (T~^(a;[2]) and 
CT~^(w[a;[l], a:[2]]) appeared after 0.031s and 5.016s respectively. The final result 

x[2][w[x[l],x[2]] - :^u;(°'^)[x[l],x[2]] - i^^(M) [^[1]^ ^[2]]) + ^^^d-O) [a;[l], a;[2]] (4.1) 
was presented immediately. 
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